home *** CD-ROM | disk | FTP | other *** search
/ MacTech 1 to 12 / MacTech-vol-1-12.toast / Source / MacTech® Magazine / Volume 05 - 1989 / 05.08 Aug 89 / Jorg's Code / mat_mul_c.c next >
Encoding:
Text File  |  1989-06-22  |  671 b   |  40 lines  |  [TEXT/MPS ]
  1. pascal void mat_mul_c 
  2.     (extended c[50][], long nc,
  3.      extended a[50][], long na,
  4.      extended b[50][], long nb,
  5.      long n1, long n2, long n3)
  6.  
  7. {
  8.     int i,j,k;
  9.     
  10.     for ( k=1 ; k <= n3; k++ )
  11.         for ( i=1 ; i <= n1 ; i++ )
  12.         {
  13.             c[i][k] = 0.0;
  14.             for ( j=1 ; j <= n2 ; j++ )
  15.                 c[i][k] = c[i][k]+a[i][j]*b[j][k];
  16.         }
  17.     
  18. }
  19.  
  20. pascal void mat_mul_c_opt
  21.     (extended c[50][], long nc,
  22.      extended a[50][], long na,
  23.      extended b[50][], long nb,
  24.      long n1, long n2, long n3)
  25.  
  26. {
  27.     register int i,j,k;
  28.     register extended sum;
  29.     
  30.     for ( k=1 ; k <= n3; k++ )
  31.         for ( i=1 ; i <= n1 ; i++ )
  32.         {
  33.             sum = 0.0;
  34.             for ( j=1 ; j <= n2 ; j++ )
  35.                 sum = sum+a[i][j]*b[j][k];
  36.             c[i][k] = sum;
  37.         }
  38.     
  39. }
  40.